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, The dynamics of a one dimensional growth model involving attachment and detachment of par- 

' tides is studied in the presence of a localized growth inhomogeneity along with anchored boundary 

, conditions. At large times, the latter enforce an equilibrium stationary regime which allows for an 

fSJ . exact calculation of roughening exponents. The stochastic evolution is related to a spin Hamiltonian 

» I ' whose spectrum gap embodies the dynamic scaling exponent of late stages. For vanishing gaps the 

CIh, interface can exhibit a slow morphological transition followed by a change of scaling regimes which 
are studied numerically. Instead, a faceting dynamics arises for gapful situations. 
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I. INTRODUCTION 



o 

1^ ' After two decades of investigations the dynamics of growing interfaces continues to be a subject of enormous interest, 
providing a framework to compare experiments, simulations and theory, let alone the wide range of applications 
encountered [1]. Despite the vast diversity of morphologies in which growing surfaces can evolve, most numerical 
, analysis and theoretical studies pointed out the onset of scaling regimes emerging at both large time and length 
scales. This enabled a classification of apparently dissimilar processes in terms of universality classes characterized 
I t \ by a set of scaling exponents which dominate the late evolution stages [1,2]. 

• A common feature associated with these processes is the possible emergence of rather slow temporal crossovers in 
^ ' which the early dynamics exhibits quite different roughening characteristics from those observed in the asymptotic 
I limit [3] . The presence of growth rate inhomogeneities or growth defects localized within small spatial regions of the 
substrate plane (columnar defects), is one of the simplest mechanisms whereby such crossovers can be observed [1,4]. 
Another possibility is realized by anchoring conditions through which nonequilibrium fluctuations are completely 
^ suppressed at the interface boundaries [5]. The main interest in those situations is in the morphological phase 
transitions that may occur at large times. In this work we investigate the change of scaling regimes accompanying 
these transitions by means of a prototype restricted solid on solid (RSOS) growth model [6,7] combining both of these 
mechanisms in one dimension. As we shall see, anchored boundaries are essential for the appearance of equilibrium 
. regimes which in turn allow for a simple calculation of roughening exponents at late stages. 
' At the phenomenological level of the Kardar-Parisi-Zhang equation [8], thoroughly studied in the continuum theory 
: of kinetic roughening, large scale morphology changes can be accounted for by assuming a growth velocity which 
■"n]" is a symmetric function of the local inclination of the interface. The analysis of this equation under a growth rate 
inhomogeneity then leads to a nonequilibrium steady state (SS) which emerges in the form of a macroscopic hill [4]. 
Here, we show that a similar scenario may also arise under equilibrium and near equilibrium conditions dominated 
by unusual scaling regimes. 

"j^ In studying the latter it is useful to consider the mean square fluctuations of the average interface height h{t) which 

yields a measure of the interface width 



w'{L,t)^jj2{[K{t)~hit)]'), (1) 

o 

O ' where the brackets denote an ensemble average over all possible evolutions of heights {/i„} forming the interface at 
time t, which grows on a substrate of size L. On general grounds it can be argued that W scales as [9] 

><■ W{L,t)^L^f{t/L^), (2) 

■ where the scaling function /(c) satisfies 

;(,)^|cC/^ for c<<l (3) 
l_ const tor 1 . 

Hence, it follows that finite systems saturate as cx X'', whereas in the thermodynamic limit the asymptotic growth 
is ruled by the exponent /3 — (/z, that is W oc . The exponent C describes the roughness dependence of the 
interface width on the typical substrate size. In turn the exponent z , often known as the dynamic exponent, gives 
the fundamental scaling between length and time. 
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In practice, Eqs. (2) and (3) yield a standard procedure which is often followed to extract and corroborate scaling 
exponents predicted by other approaches and certainly we will make use of this hypothesis as well. However, due to 
the presence of the crossovers referred to above, a complementary procedure would be needed if the former becomes 
exceptionally slow. To this aim, we will also exploit the known equivalence between the RSOS growth models already 
mentioned and a one dimensional gas of hard core particles undergoing an asymmetric exclusion process [1,6,7] 
(see Fig. 1). The idea is to cast the evolution operator of the associated master equation of this latter process 
[10] into a suitable quantum spin representation [11,12] lending itself more readily for a finite size scaling analysis. 
Since the dynamic exponent z of (2) is ultimately embodied in the gap behavior of the evolution operator (or spin 
"Hamiltonian" ) , the study of its lower spectrum can then provide information of the late evolution stages in a more 
direct manner. Evidently, this methodology along with the evaluation of the roughening exponent C - simplified 
greatly by the anchoring boundary conditions - circumvents the problem of arbitrarily long transient regimes though 
on the other hand is limited severely by the affordable system sizes. A posteriori, it will turn out that already modest 
lengths can yield clear finite size trends. This strategy will be explained in Section II and its results compared with 
those of standard techniques given in Section III. We end the paper with Section IV which contains our conclusions, 
along with some remarks on extensions of this work. 



Let us consider the dynamics of lattice aggregation models with no overhangs, including both adsorption and 
desorption of monomers at random locations of a one dimensional interface [6,7], such as that described in Fig. 1. As 
usual, on a coarse grained level of description the state of a surface at a given time is represented by a set of single 
valued functions /i„(i) measuring the surface heights at positions 1 < n < L + 1 of the growth substrate. As it was 
mentioned above, we are interested in boundary conditions that suppress completely height fluctuations at n = 1 and 
i + 1 for all times, i.e. the interface is anchored at the boundaries. For simplicity, we study the case where L is even 
and h\ = ft-L+i , whereas deposition and evaporation rates e , e' are taken uniformly throughout the system except on 
site ■§ + 1 where these probability values are respectively Cq , eg . To prevent the divergence of interface fluctuations 
in the bulk, we impose a RSOS constraint namely, — /i„| = 1 , Vn, t . Specifically, growth (evaporation) events 

hn h„ + 2 , [hn — > hn — 2], with n = 2, ... , L , occur only at local minima (maxima) of the evolving interface. These 
basic processes and their transition rates are schematized in Fig. 1. In turn, the typical configurations resulting from 
these rules at large times are displayed by the snapshots of Fig. 2. 

It is often convenient to consider the interface slope rather than the height, so hereafter we will exploit the known 
mapping between RSOS interface dynamics and quantum spin-i systems. This correspondence can be easily visualized 
in the scheme of Fig. 1. Associating the height difference s„ = hn+i — hn to an eigenvalue of the z— component, 
say, of the Pauli operator (?„ for site n, all relevant quantities of the interface, such as its width and height-height 
correlations, can be casted in terms of i- spinors. By construction, it is clear that the interface heights (relative to 
hi ) are obtained as hn = J2^=i for n = 2, ... ,L + 1 . Therefore, the anchoring condition hi = hL+i imposes the 
vanishing of the total magnetization throughout the underlying spin kinetics. 

As is well known, the probability distribution of such Markov processes is controlled by a master equation [10] 



whose transition probability rates R{s s') € {e, e', eo, ej,} , now denote the (biased) spin exchanges at which a 
generic configuration \s) = \si, ... ,sl) evolves to |s') through a single exchange of two consecutive spins. Starting 
from a given probability distribution |P(0) ) = J2s 0) |s ) , Eq. (4) can be conveniently thought of as a Schrodinger 
like representation in which the ensemble averaged state vector \P{t)) (playing the role of wave function), can be 
evaluated at subsequent times from the action of an evolution operator (or "Hamiltonian" ) on the initial state, namely 
\P{t) ) = e~^* l^'(O) ) [11]. The specific form of H can be readily obtained by introducing spin-i raising and lowering 
operators (7+, cr^ , along with spin occupation fields n = cr+cr" . It is then straightforward to show that the stochastic 
dynamics of Eq. (4) is accounted by the operator 



II. SPIN REPRESENTATION 




(4) 
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where e„, are shorthands denoting respectively 

{eo, e'o for n = ^ , 
(6) 
e, e' otherwise , 

while the anchoring condition hi = Hl+i confines the dynamics within the subspace ri„ = L/2 . The biased hop- 
ping terms of (5) clearly describe the original growth-desorption events (see Fig. 1), while conservation of probability 
requires in turn the appearance of the remaining (diagonal) particle-vacancy correlators. We address the reader to 
Ref. [12] for a more detailed derivation in related systems. 



A. Detailed Balance 



Before continuing with an alternative spin representation more suitable to study dynamical aspects at large times, 
we pause and consider the SS of Eq. (5) along with its implications on the interface character. 

Given two spin configurations \s) = |... , s„, — s„, ... ) , |.s' ) = |... , — ,s„, s„, ... ) , differing at most in the state of two 
neighboring n, n + 1 locations, evidently detailed balance probabilities in (4) will hold provided that 



P(s)4 = Pis') en, if {s'\H\s) = -e'„, 
P(s)e„ = P(s')4, if {s'\H\s) = -en. (7) 

This can be readily satisfied defining a hard-core particle (up spin) potential 

V{n) = E In j = n In (1) + In (^^) (n - L/2) , (8) 
through which the equilibrium distribution is simply obtained as 



P(si, ... , Si) oc exp- 



^J2vin){l + sn) 



(9) 



When e — e', these probabilities further enable us to construct the partition function (normalization constant), height 
profiles (spin densities), as well as the spin correlation functions needed to derive the equilibrium interface width. 
For e ^ e' a rather involved recursive relation in the particle number can be obtained for all these quantities, but 
its analytic solution is not reachable by standard means [13]. However, this case does not yield a rough interface as 
statistical fluctuations become exponentially suppressed in time (see Section III B). 

If e = e', the step function potential (8) permits to divide the system into two independent regions [1, ^], [^ -|- 1, L] 

with (^) configurations having ^ — m and m-particles respectively (0 < m < ^). Hence, the partition function 
normalizing the above SS distribution is given by 



m=0 



(10) 



where r = Eq/eo • Using analogous arguments, we can also obtain the reduced partition function Zi which arises from 
the occupation of a given site i, and evaluate the spin density (a|) = 2 Zi/Z — 1 . This results in 

for i<f , 

Zi = { T_\ (11) 
E LtJ(*-^)r-+\ otherwise, 

m=0 

which implies an equilibrium shock profile stemming entirely from the inhomogeneous potential (8) at finite particle 
densities. In particular, for L — > oo the analysis of Eqs. (10) and (11) yields the following discontinuity 



«) = [l-2e(n-L/2)](s), {s) = ^^^-^. (12) 

V eo + V^o 
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Thus, in the hight representation, so long as eg/eo < 1 (> 1 ) Eq. (12) entails the sideways growth of a hill (valley) 
whose sides at large times are tilted by an amount of ±(.s) relative to the substrate. 

To determine whether this morphology is actually rough, we focus attention on the equilibrium height fluc- 
tuations 5n = {h'^n) ~ (^n)^ • Consequently, first we evaluate the spin-spin correlations involved in (/i^) = 
(n — 1) + 2 ^ {af aj) . Thus, once more we recur to the combinatorial reasoning and calculate the reduced partition 

i<j<n—l 

functions Zij resulting from the occupation of two specific sites i,j . After some elementary steps we obtain 



Zij — { 



i-2 

E (*-')(*) for i<i<^, 

m=0 

*E (*;})( t^)r-+\ for i<f <i, 

m=0 

ijjit')^"'^'^ othrewise, 

m=0 



(13) 



from which the required spin correlations are computed as (cr| it|) = 4:Zij/Z — ( (erf) -|- (cr|) + 1 ) . In the large size 
limit the analysis of Eqs. (10), (11) and (13) ultimately yields a rough interface (see leftmost snapshot of Fig. 2), 
whose height fluctuations (canceled at the boundaries), result distributed as (5„ oc (n-1) ( 1 - 1^=1 ) . This simple form 
contrast with that observed in nonequilibrium SS of anchored self organized interfaces in which fluctuations in the 
upper part of the hill are substantially reduced [5] . 

The above correlations can also characterize the saturation width referred to in Section I. Specifically, in the spin 
representation it can be easily checked that W"^ may be rewritten as 



L2 



6L 



+ 



2 

Z2 



z{L-j){ata]). 



In Fig. 3 we display the size dependence of W for several values of fo/eo. 
bears the discontinuous character of Eqs. (11) and (13). More precisely, 



(14) 



It turns out that the roughness exponent C 



1 , if 7^ eo , 
1/2 , if e'„ = eo. 



(15) 



Often, a value of = 1 is special because it signals that the assumption of a well defined average orientation of the 
interface (parallel to the substrate plane), becomes inconsistent. Certainly, this is in line with the tilt obtained in Eq. 
(12). For c'q = eo the conventional (diffusive) roughening is recovered; here the tilt vanishes and the orientational 
fluctuations at large scales estimated, for example, as W{L,t oo)/L, decrease with L. We shall revisit this point 
later on in Section III A. 



B. Self Adjoint Representation 

As is known [10], detailed balance guarantees the existence of a representation in which the evolution operator 
(5) is self-adjoint. Although an exact solution of the (real) i7-spectrum in the thermodynamic limit seems unlikely 
irrespective of its representation, at least a self-adjoint description can facilitate the numerical analysis of a finite size 
scaling approach. Specifically, one can readily find a similarity transformation to map (5) into an hermitian matrix, 
and thereafter obtain the lower eigenmodes dominating the asymptotic kinetics via recursion-type algorithms, e.g. 
the Lanczos technique [14], appropriate to study fair system sizes. 

To this aim and with the aid of the particle potential introduced in (8), we rotate the above operator around the z 
spin direction using a pure imaginary site dependent argument ip{n) 



This rotation is produced by the non-unitary similarity transformation U = e 
turn results in the direct product 



-iS 



e 4 



V{n) 



e 4 







(16) 

with S = \ Ylin ^(^) ^ra ' which in 



(17) 
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While the diagonal terms of (5) remain unaltered by U, it is straightforward to show that 

C/„a±C/-i = eT'^(")a±. (18) 

From this latter transformation, one can immediately verify that the rotated (self- adjoint) operator H = U HU~^ 
becomes an open XX Z ferromagnet with a defect coupling under local magnetic fields, namely 

1 

W = -2 E Jn + + A„«<+i - 1)] 



h{c7l - cjl) - {ho-h){al - al,^). (19) 



where 



A„ = (e„ + 4)/v/4e„ 

/io = (eo - eo)/4, (20) 
h = {e- e')/4, 

with e„, ejj taken as in Eq. (6). Thus, we are left with a diagonalization problem which, to some extent, is now 
controllable by standard recursive techniques (Section III). 

For the sake of completeness it is worth pointing out that the similarity transformation (17) also enables us to 
obtain the SS distribution (9). In fact, exploiting that is a stochastic operator, we can express its left SS as 
an equally weighted sum of all accessible configurations [10], i.e. ^\ = (0| X^„^...„^ cr~^ ■■■^nL ' where (0| denotes the 

ferromagnetic down spin state. Hence, by construction |'^o) = U~^\'ip) is a (unnormalized) ground state of H , and 
therefore the SS distribution in the initial if-representation is constructed as 

£/-|^)« Y: e-^(-)...e-^K^^<...a+jO), (21) 

ni...nL ^ 

2 

thus recovering the equilibrium distribution (9). 

Returning to the dynamics, the average value of a diagonal quantity VV - such as the "width operator" involved in 
Eq. (14) - varies according to {tp\ VVe"^* l-P(O)) [12]. Since VV is invariant under U, it is a simple matter to check 
that in the self-adjoint representation (W) reads 

(W)(L,i) = (^o|W|Vo) + e-^*(^o|W|^A)(^A|C/|P(0)), (22) 

where {|V-'a)} is a complete orthonormal set of eigcnstatcs of H (all with A > 0). As expected, the role of initial 
conditions becomes irrelevant near the equilibrium regime. If the spectrum gap vanishes in the thermodynamic limit, 
the width approach to equilibrium will involve arbitrarily large times for sufficiently large systems. In those situations, 
finite size scaling analyses of the first excited levels A^, would then provide the dynamic 2;-exponent ruling over the 
late roughening stages referred to in Section I. 

A distinctive feature arises when all components of the total angular momentum S = ^ a„ are preserved by H, 
namely for e = e' and = • Since W just involves operators of the form cr| [see Eq. (14) ], then rather restrictive 
selection rules hold for its matrix elements in (22). Specifically, given that {|V'a)} can be classified according to the 
total spin 5 [ = S'(S'-|- 1) ], the non- vanishing contributions to Eq. (22) come only from states having 5=^ — 1 
and ■§ — 2 [15]. What should be emphasized here is that as soon as [H, S] 7^ the effective density of states, partly 
responsible for the temporal asymptotic behavior of (22) when L ^ 00 , is drastically modified as these selection rules 
no longer apply. We will come back to this issue within the numerical context of Section III A. 



III. NUMERICAL RESULTS 



To explore the dynamical consequences of these arguments we have carried out Monte Carlo simulations as well as 
finite size scaling analyses of the RSOS model referred to above, for a variety of situations. First we focus attention 
on the subcase e = e' where the roles of eo and Cq are clearly interchangeable, so we restrict the analysis to, say, 
Eo/eo < 1 • In this situation the interface actually roughens and exhibits two different scaling regimes accompanying a 
large scale morphological transition. The discussion of different bulk probability rates is addressed at the second part 
of the Section. It will turn out there that fluctuations decay very rapidly and prevent the roughening of anchored 
interfaces at large times. Instead, a faceting dynamics will emerge regardless of the inhomogeneity growth rates 
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A. £ = €' 



Starting from an initially flat configuration we studied the evolution of the interface width and monitored the height 
profiles at different growth stages. Fig. 4 displays the width behavior obtained for eQ/eo = 0.5 and 0.8. The rather slow 
crossover (particularly for 0.8, where it is only incipient), deterred us of using larger substrates, though preliminary 
simulations averaged over few evolution samples indicated similar trends. The early growing stages support a power 
law growth W oc extended over more than four decades with an exponent /? ~ 1/4 . This typical diffusive behavior 
is accompanied initially by a height profile which is almost parallel to the substrate, except in the neighborhood of 
the growth inhomogeneity. As is shown in the inset of Fig. 4, on approaching the asymptopia however, the slopes of 
the hillsides steepen until they reach the equilibrium tilt alluded to in Eq. (12). This progressive orientation departure 
signals a large scale morphological transition which in turn is also reflected in the increase of the growth exponent. 
In the late dynamic stages, this can be well fitted by a value of /? ~ 1/2 for nearly two decades. 

An alternative determination of this rather peculiar value [16], can be implemented by resorting to the phenomeno- 
logical scaling assumption referred to in Eq. (2) along with the findings and arguments of Section II. Notice that this 
provides two independent numerical procedures to estimate the dynamic z-cxponcnt needed for the knowledge of (3. 
On one hand, the former can be calculated by studying the finite size behavior of the first excited levels of H (in 
principle, just the lowest will do), so we diagonalized it exactly via a recursion type Lanczos algorithm [14] applied 
on the zero magnetization subspace. The huge dimensionality of this sector, growing as (j^yg)' 'wd^ as the lack 
of translational symmetry of the evolution operator, limited our computations to chain lengths of up to 24 spins. 
Nevertheless, they proved to be sufficient for a fair estimation of the spectrum gap. We direct the reader to Fig. 5 
which suggests a decrease cx for the gap and other excited levels consistent with a common value of z = 2 . On 
the other hand, an independent evaluation of this exponent can be attained by studying the scaling behavior of the 
interface width [Eq. (2) ]. Fig. 6 exhibits the results of our simulations for growth substrates of 2500 and 5000 heights 
on approaching their saturation regimes. Here, the data collapse was obtained upon setting (2;, C) ~ (2, 1) , which 
confirms not only the Lanczos estimation but also corroborates the roughness exponent quoted in Eq. (15). Thus, 
from Eq. (2) it follows that the fast roughening behavior already observed in Fig. 4, is now recovered by the ratio 
C/z. 

Since the Lanczos analysis continues to yield values of z = 2 holding up to the homogeneous situation (as it should), 
a natural question one can pose is therefore: through which feature does the dynamics render a completely different 
roughening behavior at large times as soon as cq c'q 7 In an attempt to provide a plausible explanation for the 
appearance of this abrupt change [see also Eq. (15) ], we resort to the observations given by the end of Section II B. 
The; inset of Fig. 5 displays the lower part of the 7i-spectrum for both inhomogeneous and homogeneous situations. 
As stated above, the latter involves at most L{L — l)/2 contributing levels with S = L/2 — 1, L/2 — 2 , while in the 
former case the sum of Eq. (22) becomes much denser as [H, S ] ^ and an exponential number of new states arises. 
Of course, the new matrix elements might eventually change from zero in a continuous manner, but the density of 
states (a measure of which is given by the inverse of the levels spacing), varies abruptly. Moreover, some of the low 
lying c;xc;itations c;ontrolling the asymptotic regime (t — > 00 holding t/L^ ^1), suggest a rather narrowly peaked 
structure which is entirely absent for cq = c'q . 

Also, it is interesting to examine whether the large scale morphology transition embodied in the slow temporal 
crossover of W affects the usual scaling hypothesis of Eqs. (2) and (3). Thus, we turn to the early dynamic scaling 
of W shown by the inset of Fig. 6. In contrast to the super-diffusive growth observed at late stages, here the 
data collapse arises by setting standard diffusive exponents {z, ~ (2, 1/2) which in turn yield a scaling function 
cx (t/L'^y/^ , in agreement with the early /? exponent measured in Fig. 4. Hence, combining the late and early scaling 
regimes it follows that 

WiL,t,r) = L^Mt/L^), (23) 

where r is a crossover time which depends solely on Cq/cq (eventually diverging in the limit Cg ^ eo ), and /r(c) is a 
universal function defined over three different scales as 

fc^/^ for c«;r/L2, (C = 1/2) , 

/r (c) ~ <^ c 1/2 for r/L2 « c « 1 , (C = 1) , (24) 
I const for c ^ 1 . 

In what follows we address finally to the e ^ e' situation. 
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B. ey^e' 



A quick glance at the evolution of the interface width displayed in Fig. 7 might render the (wrong) impression that 

the general case e 7^ e' , eg 7^ bears similar characteristics. However this time the growth exponent exhibits an 
asymptotic value of /3 ~ 3/2, which is not understandable in terms of conventional kinetic roughening theories [16]. 
In fact, as we shall see below, the interface does not roughen. 

Let us first provide a simple explanation for this large value of (3 by means of the following heuristic considerations. 
For clarity of argument, assume vanishing desorption rates and e = eg. Starting from a flat configuration, say with 
h2k = 1, h^k-i = , a deterministic dynamics arises as (L — 2)/2 deposition attempts occur on the initial (L — 2)/2 
interface minima. Next, we are left with h2k = 1, h2k-i = 2, hi = hz+i = and {L — 4)/2 contiguous minima 
over which new (L — 4)/2 depositions will be once again deterministic. By iterating this argumentation t-times, 
this dynamics leads to a configuration resembling a truncated pyramid (see central snapshot of Fig. 2, though for 
for a non-deterministic situation, i.e. e' 0). Specifically, there will be {L — 2t)/2 contiguous minima in between 
ht+2 =t+l and hL-t-i = ht+2 ■ Using the common definition of Eq. (1), it is then easy to verify that the "width" 
of such configuration is simply 



W{L,t) 



0{t/L)] , i.e. 13 = 3/2. 



(25) 



The above argument describes rather a faceting process (terminating at f ~ L/2 ), which strictly applies for e' = . 
Certainly, as soon as e' > the dynamics is no longer deterministic, no matter how small e' is. However, for e' <C 1 
the early desorption attempts become gradually unsuccessful as the active region of the interface, i.e. the number of 
available minima, decreases inasmuch the sideways region is increasingly jammed (Fig. 2). Thus, at large times a 
process similar to a faceting dynamics might be expected, at least for small bulk desorption rates. In fact, for e > e' 
and eo > o^^'" numerical simulations confirm these considerations. No differences were observed between sampling 
histories at large times, so fluctiiations become asymptotically negligible. For eq < situation is similar though 

another faceting process shows up around the inhomogeneity, as displayed by the height profiles in the inset of Fig. 7 
(see also rightmost snapshot of Fig. 2). Of course, for e < e' the roles of eo and Cq are interchanged. Ultimately, the 
whole process approaches a non-fluctuating pile of slope ±1, so long as e' ^ e . 

To provide an alternative understanding of this fast fluctuation decay for generic rate values, we recur once more 
to the analysis of the spectrum of the evolution operator (19). The gap and levels obtained for the sizes within our 
reach are shown in Fig 8, but in contrast to the e = e' situation, here the finite size trend of these quantities needs 
further analysis. To this end, we studied the 7i-spectrum within the subspace = L/2 — I (single spin excitation), 
which corresponds to the much simpler anchoring case — hi = L — 2 . For this sector, it is straightforward to 

check that Eq. (19) reduces to the tridiagonal matrix 



H = 



e' 7 
J e + e' 7 







7 £ + 4 7o 
7„ e'-heo 7 



7 e + e' 



(26) 



where 



7o 



(27) 



Evidently, this just constitutes the simplest approximation to the many body problem of = . Nevertheless, the 
comparisons of Fig. 8 indicate that the eigenvalues of (26) yet provide an excellent estimation of the actual gap and 
other excited levels obtained through the Lanczos scheme. (It is worth pointing out in passing that an excellent fit of 
these quantities was also found for e = e' ). Using the spectrum gap of the single spin approximation, the inset of Fig. 
8 strongly suggests that the same gap will persist for Ai in the thermodynamic limit of (19). Similar gapful results 
were obtained for other values of Cq/eo • Thus, fluctuations would be suppressed at large times, which is in line with 
the almost invariant values of W observed over many sample histories. 
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Also, at the gap level the density of states diverges as L in the simplified version of the problem. For = 
however, one might eonjecture a much stronger divergence, probably growing like ^ , as the number of levels 
between two single excitations tends to increase exponentially with the system size (at least for the small lengths at 
hand). This would leave us with a saturation time oc i in Eq. (22), which on the other hand would be in agreement 
with the termination time of the faceting process idealized above. In fact, the numerical simulations displayed in 
Fig. 9 lend further support to these speculations. Clearly, these results exhibit both a saturation time oc L , as well 
as an asymptotic scaling regime consistent with Eq. (25), i.e. W/L oc {t/L)^^^ , during which both height and width 
fluctuations are absent. 

In contrast, at early stages the interface displays typical roughening features. Specifically, the inset of Fig. 9 
exhibits a diffusive scaling regime W/VL = f(t/ L^) ,'m. turn corroborated by the growth exponent /? ~ 1/4 obtained 
in larger systems (Fig. 7). This strong departure from the faceting description occurs on temporal scales smaller than 
a crossover time which turns out to decrease when e'/e — > (irrespective of Cg/eo), but eventually diverging in the 
limit e' ^ e , eg ^ Co . 

IV. CONCLUSIONS 

We have analyzed the characteristics of both early and asymptotic dynamics of one dimensional anchored inter- 
faces under a growth inhomogeneity. There are two sets of results related respectively to equal or different growth- 
evaporation rates in the bulk. 

For e = e', even the slightest departure from the homogeneous eo = c'q situation is able to produce finite interface 
tilts as well as huge temporal crossovers. The problems posed by the latter have been bypassed studying separately 
the roughness C, and dynamic z exponents, the evaluation of which was significantly simplified by the appearance 
of equilibrium SS. In analyzing finite size scaling trends of the spectrum gap via the Lanczos method, we found no 
changes with respect to the homogeneous situation, i.e. z = 2 . Consequently, it was argued that the breaking of 
the full rotational invariance of the evolution operator is iiltimately responsible for the emergence of a much heavier 
density of states accounting both for the discontinuity of the roughness exponent C, (or alternatively, for a different 
asymptotic growth exponent /3), as well as for the rise of a new scaling regime at large times. These expectations 
were confronted independently with standard numerical simulations monitoring the evolution of the interface profile 
and width. At early stages, the latter exhibits a diffusive scaling regime having basically a non-tilted profile (except in 
the inhomogeneity neighborhood), but progressively approaching a final regime consistent with our scaling exponents 
C, = 1, z = 2 . However for eo ^ eg , in practice this new regime might occur at a time so large as to render it 
numerically unobservable. 

For e 7^ e' the situation is entirely different. Here, the spectrum gap docs not vanish in the thermodynamic limit 
regardless of the inhomogeneity rates, and fluctuations between evolution histories at large times become negligible. 
This confirms an heuristic description (in turn, tested independently by simulations), suggesting that the asymptotic 
dynamics becomes almost deterministic. We may also think of a synchronous discrete time process in which a 
randomly chosen finite fraction, or possibly all of the growth sites, are simultaneously updated in a single time step. 
One characteristic feature of such synchronous models is the occurrence of faceting transitions at large times [1], which 
also turned out to be the case here. In contrast, at early stages the interface actually roughens following a typical 
diffusive pattern accompanied by a standard scaling regime. 

The analysis of nonequilibrium asymptotic situations, even for d = 1 , might become rather involved. In this sense, 
it will be interesting to elucidate whether a direct evaluation of C, could be achieved using the matrix approach to 
the asymmetric exclusion process [17] with both injection and ejection of particles at the boundaries, including one 
or more hoping defects (that is, unanchored boundaries and growth inhomogeneities in the height representation). 

Higher dimensional extensions of this study would be clearly desirable and more realistic. However, the analysis of 
the corresponding quantum spin analogy should involve a projector operator to discard all those spin configurations 
having magnetic loops, i.e. cr^ 7^ , with r on a given closed path. Otherwise, the mapping would no longer 
represent an interface. The issue as to whether or not such ideas are actually practical in d > 1 remains quite open. 

The author acknowledges support of CONICET, Argentina. 
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FIG. 1. Schematic representation of monomer deposi- 
tion-evaporation onto a RSOS interface with anchored bound- 
aries at hi = hh+i The equivalent spiu-i (sn = /tn+i — hn), 
or hard core particle dynamics involves a left (right) particle 
hopping with rate e (e') for monomer adsorption (desorption). 
The corresponding rates for the inhomogeneity at hhj^^ are 
€o and ej). 




FIG. 2. Possible evolution scenarios. Typical snapshots for 
L = 10^ using e'/e = 1 with e'o/eo = 0.5 after t = 10^ steps 
per height (left); e'/e = e'o/eo = 0.5 at t = 500 (center); and 
e'/e = 0.5 with e'o/eo = 5 at t = 500 (right). For e ^ e' , 
fluctuations are progressively reduced on their way to the pile 
configuration denoted by dotted lines with slopes ±1 . 
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FIG. 3. Equilibrium interface width for e/e' = 1. Up- 
per and lower dashed curves denote respectively the cases 
^o/^o = and 1, whereas solid lines going downwards stand 
for eo/eo = 0.1, 0.3, 0.5 and 0.7. 
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FIG. 4. Growth of interface width for e/e' = 1 using 
€o/co = 0.5 (circles) and eo/eo = 0.8 (squares showing an 
incipient asymptotic deviation) for L = 10'' averaged over 
200 histories. The early and late slopes of dashed lines, are 
respectively 3 ~ 1/4 and 1/2. The inset displays the height 
profile evolution for eo/eo = 0.5 averaged over 2000 histo- 
ries with L = 1000, at t = 1.5 x 10*, 6 x 10*, 2 x 10^ and 10^ 
[dotted line at the top following closely the tilt s of Eq. (12) ]. 
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FIG. 5. Finite size behavior of three lowest excited levels 
of the evolution operator (5) for e/e' = 1 and eo/eo = 0.5. 
Solid lines have slopes z = 2 . The inset compares the lower 
spectrum of this case (L = 20 , 10 spin-excitations, framed 
at the left), with the effective levels of the regular situation 
(eo/ef) = 1, L = 20, 2 spin-excitations). 
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FIG. 6. Asymptotic finite size scaling regime of the in- 
terface width averaged over 200 histories using e/e' = 1 , 
e'o/eo = 0.5, for L = 5000 (squares), and L = 2500 (cir- 
cles). The saturation values (horizontal line), coincide with 
those of Eqs. (13) and (14). The dashed line is fitted with the 
slope calculated from the data of Figs. 3 and 5 (i.e. I3 = (/z). 
The inset exhibits an early scaling regime which behaves dif- 
fusively (/? = 1/4), alike the initial data of Fig. 4 
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FIG. 7. Evolution of the interface width for e' /e = 0.5 , 

^o/^o ~ 5 and L = 10® . The initial and final slopes of 
dashed lines are /3 = 1/4 and 3/2 respectively. From top 
to bottom the inset exhibits the profile of 10^ heights at 
t = 5 X 10", 2 X 10", 10* and 4 x 10^ . Both width and height 
fluctuations become negligible at large times. At the width 
level, results of different /eo values closely follow each other 
in all evolution stages. In turn, for eg < eo the early profile 
has no tilt around the inhomogeneity. 
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FIG. 9. Late scaling regime of W for e'/e — eo/eo ~ 0.5 
using L = 10'' (circles) and L = 2 x 10** (squares). The 
dashed line slope and scaling form are both consistent with 
the faceting process conjectured in (25), whereas differences 
between sample histories become gradually negligible. In 
contrast, the inset results, averaged over 200 histories, in- 
dicate an early scaling regime which is typically diffusive 
(C = l/2, « = 2, /? = C/a:). 
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FIG. 8. Finite size behavior of lower excited levels in the 
single spin approximation (5^ = L/2 — 1 ) of operator (26) for 
e'/e = eo/eo = 1/2 (solid lines). Circles and triangles stand 
respectively for levels Ai and A2 of Hamiltonian (5) with 
= . Above A2 , further collective excitations (not shown) 
appear between successive solid lines. The inset suggests a 
power law convergence of Ai(L) (circles, =0), towards 
the gap g obtained in the main panel. 
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